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In this work the dynamics of a freely jointed random chain with small masses 
attached to the joints is studied from a microscopic point of view. The chain is 
treated using a stringy approach, in which a statistical sum is performed over all two 
dimensional trajectories spanned by the chain during its fluctuations. In the limit 
in which the chain becomes a continuous curve, the probability function for such a 
system coincides with the partition function of a generalized nonlinear sigma model. 
The cases of open or closed chains in two and three dimensions are discussed. In 
three dimensions it is possible also to introduce some rigidity at the joints, allowing 
the segments of the chain to take only particular angles with respect to a given 
direction. 



I. INTRODUCTION 

In this paper the dynamics of a random chain subjected to thermal fluctuations at fixed 
temperature T is discussed. This problem is usually treated phenomenologically, regarding 
the fluctuations of the chain as a stochastic process which may be described with the help of 
Langevin equations or, alternatively, of Fokker-Planck equations (if . This approach leads to 
the well known models of Rouse Q and Zimm ^ which allow a satisfactory understanding 
of the main properties of polymers in solutions. One major drawback of these coarse grained 
models is, that they suffer from the presence of rigid constraints. The Rouse and Rouse- 
Zimm equation consider only chains bead spring models of chains, where the local spring is 
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infinitely extensible. In the framework of continuous models the Rouse equation for example 
is nothing but the stochastic equation (Langevin equation) for the classical Wiener measure, 
which yields paths which do not have a well defined tangent vectors These problems 
have been tackled by various attempts, see e.g. However, the correct use of rigid 

constraints in (stochastic) dynamics requires some mathematical effort 0,(3], in contrast to 
the static cases where rigid constraints can be implemented by Dirac delta functions in the 
partition function. 

However, it seems to be too difficult to extend these approaches, therefore we follow here 
some different routes. Moreover the use of dynamic constraints for (stochastic dynamics 
equation) is always problematic. 

Here the dynamics of random chains is considered from a microscopic point of view. The 
chain is represented as a set of freely jointed segments of fixed lengths. Small masses are 
attached at the joints. The goal is to construct the probability function \l'[Confj, Confj] of 
this system. Roughly speaking this function measures the probability that a chain starting 
from a given initial conformation Confj at the time ti arrives to a given final conformation 
Confj at time tf. During its motion from the initial to the final conformation the chain 
spans a two dimensional surface. In this way one obtains a stringy formulation of the chain 
dynamics, which is based on path integrals 8]. The chain "world-sheet" is a rectangle whose 
sides are given by the length of the chain and the time interval At = tf — t^. 

The basic ingredient in the construction of the probability function \E'[Conf/, Confj] is the 
energy of a discrete freely jointed chain with segments. While it is easy to add external 
and internal interactions acting on the segments of the chain, the complicated form of the 
kinetic energy poses serious obstacles to the possibility of computing the probability function 
analytically. The kinetic energy is in fact nonlinear and has an intrinsically non-Markoffian 
nature. This is due to the fact that the velocity of each segment in the chain depends on 
the velocities of the other segments. In this situation, it is difficult to isolate in the energy 
a linear term which could give raise to a propagator, so that even perturbative calculations 
are not allowed. Luckily, the expression of the kinetic energy simplifies considerably after 
performing the limit in which the chain becomes a continuous curve. The final form of the 
probability function which one obtains in this limit closely resembles the partition function 
of a nonlinear sigma model Q]. The difference with respect to the standard nonlinear sigma 
model is that in the latter case the modulus of the fields is constrained, while in the present 
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case the constraint involves the modulus of the derivative of the fields with respect to the 
arc-length of the chain. This constraint is related to our assumption that the lengths of the 
segments is fixed. 

Within our formalism it is possible to add constraints to the trajectory of the chain 
with the help of Dirac 5— functions as in the case of the statistical mechanics of random 
..ains 0. He. we have ccs.de.ed ,„st the simplest exa,„p.e of co„st.a,„t. na.ely the 
requirement that the trajectory of the chain is closed. Besides, it is easy to treat non- 
homogeneous chains, in which both lengths of the segments and the masses located at the 
points are arbitrary. The probability functions is constructed for two and three dimensional 
random chains. In three dimensions one may introduce further constraints, which fix for 
instance the lengths of the projections of the segments on the z axis. In this way we are able 
to discuss also rigid chains, in which the segments are allowed to form only a given angle 
with respect to the z axis. 

The material presented in this paper is divided as follows. In Section II the expression of 
the kinetic energy of a discrete chain with segments is derived in two dimensions using 
a recursive method. The rules of the passage to the continuous limit are established. After 
this limit is performed, important simplifications occur. To further simplify the problem, 
one of the ends of the chain is fixed at a given point. Moreover, it is supposed that the 
distributions of masses and segment lengths along the chain are uniform. It is shown that 
under these assumptions the kinetic energy of the continuous chain may be written in terms 
of free complex scalar fields subjected to a constraint. The origin of this constraint is 
the requirement that the lengths of the segments are constant. In Section III the classical 
solutions of a free chain are studied. In Section IV a stringy approach to the dynamics of the 
chain in two dimensions is established. The probability function is constructed using path 
integrals. The resulting model closely resembles a nonlinear supersymmetric model. The 
problem of fixing various boundary conditions, including the case of closed chain trajectories, 
is discussed. A perturbative approach based on the expansion of the fields describing the 
statistical fluctuations over a classical background is provided. In Section V the results of 
Sections II-IV are extended to three dimensional chains. In Section VI chains with constant 
bending angles are investigated. Finally, in Section VII our conclusions are presented. 
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II. THE ENERGY OF A FREE CHAIN 

Let us consider a chain of segments of fixed lengths I2, ■ ■ ■ ,In in the two dimensional 
plane. Each segment Pi+iPi is completely specified by the positions of its end points Pj+i 
and Pi. In cartesian coordinates (x, y) these positions are given by the radius vectors: 

rj = {xi,yi) i = l,. ..,N (1) 

The segments are joined together at the points Pi, where 2 < / < — 1, see Fig. ^ while 
Pi and Pn are the ends of the chain. Moreover, at each point Pi, with z = 1, . . . , A^, a mass 
rrii is attached. In the following we restrict ourselves to the case of a free chain. We will see 
below that the addition of interactions is straightforward. 

To compute the kinetic energy of the above system, it is convenient to pass to polar 
coordinates k, Lpi as follows: 

n n 

Xn = ^liCOSLpi yn = ^lism(pi {n=l,...,N) (2) 

i=l 1=1 

The (fii is the angle formed by segment i with the y— axis, see Fig.^ All the radial coordinates 
li are constants for i = 2, . . . , N. The only exception is the length li{t) which denotes the 
distance of the point xi,yi from the origin. Since this distance is not fixed, li = li{t) is 
allowed to vary with the time t. From Eq. (0) the velocity components of the n— th segment 
may be written as follows: 

n-l 

in = -"^h^PiSiinpi- ln(Pn sin tpri + k cos ipi {n = 2,...,N) (3) 

i=l 
n-l 

i/n = '^h'Pi COS ipi + ln<^n COS Lpn + k sin Lpi {u = 2, . . . , N) (4) 

1=1 

Xi = —li^pi sin Lpi + li COS (fi (5) 
iji = cos + /i sin (6) 

Separating the contribution coming from the first n — 1 variables, the kinetic energy X„ of 
the n— th segment can be expressed in terms of the kinetic energy Kn-i of the (n — 1)— th 
segment: 

Tn Tfl n—l 

Kn = —Kn-1 + -;^ll'Pn + "^n XI ^nk'Pn'^i COs((/7i - + mJn'^Jl sin((/7i - (7) 

1^n—l ^ i—i 
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FIG. 1: A chain with TV segments. Let us note that the end point Pi is not bound to stay at a fixed distance 
with respect to the origin of the cartesian reference system. 



It is possible to solve the above recursion relation to find an expression of Kn- If we do that, 
at the end the total kinetic energy of the discrete chain: 

N 



becomes: 



Ty-2d 

disc 



N n-1 



(8) 



n=l 



disc ^ ^ (^iV^l + ^l) + h^l COs((/?„_fc+i - ifi) 

^ n=l k=l 

N n-1 

+ ^Jn-k+l^n-k+l Sin(v9i - ipn-k+l) 

n=l k=l 

N n-1 ^ N n-ln-k 

+ YY^l-k+l^'fl-k+l + YYY^nla-k+lk^n-k+mCOsiipn-k+l- ^t) (9) 
n=l k=l n=l k=l i=2 



where M = Yln=i "^n is the total mass of the chain \ll 

We wish now to perform the limit in which the chain of segments becomes a continuous 
system jl^. To this purpose, it is convenient to consider the indices i,k,n,... appearing 
in Eq. Q as discrete variables taking values in a one dimensional lattice with points. 
Quantities carrying the index i may be interpreted as functions of i. Their variations 
Afi are given by: A/j = /j+i — /j. Clearly, AjZ = 1, i. e. the spacing between two 
neighboring points in the lattice is 1. In order to proceed, we rescale the distances in the 
lattice in such a way that the spacing in the new lattice will be a. To this purpose, we 
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perform the transformations i — > Sj, /j — > fi^i) where the new variable has variation 
Asi = Sj+i — Si = a. The next step is to compute the kinetic energy of Eq. Q in the hmit 
— y oo, a — > 0, while the product Na remains finite, let's say Na = L, where L denotes 
the total length of the chain. Clearly, in this limit the right hand side of Eq. will diverge 
unless we suppose that the masses rrii and the lengths k of the segments are going to zero 
in a suitable way. Reasonable assumptions are: 

k — > l{si) = api{si) iTLi — > m{si) = apm{si) (10) 

where pi{si) and Pm{si) are respectively the distribution of length and of mass within the 
chain. Here we allow for segments of different lengths l{si) and for points Pi of different 
masses To be consistent with our settings, the distributions pi{si) and Pm{si) must 

be normalized as follows: 

N N 

J2 Piisi)Asi = L Pm{si)Asi = M (11) 

i=l i=l 

At this point we are ready to pass to the continuous limit. Functions of discrete variables 
will be substituted with functions of continuous variables, while sums will be replaced with 
integrals according to the following rules: 

f{s,)-^fis) 5:As, ^ / ds (12) 

i=l -^0 

After a few calculations one finds: 

K^c — K^' (13) 

where 

+ (pi{t)li{t) / dspm{s) / dupi{s — u)(p{t, s — u) cos{(p{t, s — u) — ipi{t)) 
Jo Jo 

+ h{t) dspm{s) dupi{s - u)if'{t,s - u)sm{ipi{t) - ip{t,s - u)) 
Jo Jo 

pL PS rs—u 

+ dspm{s) / dupi{s ~ u) I dvpi{v)(p{t, s — u)(p{t,v) cos{(p{t, s — u) — (p{t,v)\14:) 



and 



dspiis) = L ["^ dspmis) = M (15) 

Jo Jo 



Eq. ()14|) may be simplified by performing in the integrals in du the following change of 
variables: 

u' = s — u du' = —du (16) 

and then using the formula: 

ds r du'fiu') = ds(L - s)f(s) (17) 
JO Jo Jo 

which is valid for any integrable function f{s). As a result, we obtain: 

+ (pi{t)li{t) I ds{L - s)pm{s)pi{s)(p{t,s)cos{(p{t,s) - (pi{t)) 
Jo 

+ h{t) ds{L - s)pm{s)pi{s)(p{t,s)sm{(pi{t) - ip{t,s)) 





+ / ds{L - s)pm{s)pi{s) [ dupi{u)ip{t,s)ip{t,u)cos{(p{t,s) - ip{t,u)) (18) 
JO Jo 

For simplicity, we suppose the length and mass distributions in the chains are uniform. As 
a consequence, we put: 

M 

Plis) = 1 pm{s) = — (19) 

The first of Eqs. (fT^ implies (see Eq. (fTUj) ): 

li = a i = 2,...,N (20) 

Remembering the definition of the length distribution in Eq. (fTUI) . it is easy to realize that 
the first of Eqs. p9|) implies that all segments of the chain have the same length. As a 
further simplification, we will study the case in which the point Pi is fixed, so that 

h = ipi = (21) 

From the above assumptions and from Eq. (fTHj) . we find that the total energy Hq of the 
ideal chain is given by: 

Hq{(p) = — / ds{L — s) / du(p{t, s)(p{t,u) cos{(p{t, s) — (p{t,u)) (22) 
L Jo Jo 

In view of a future path integral treatment of the dynamics of the chain, it will be convenient 
to use cartesian coordinates, because of the complications involved with polar coordinates 
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in the path integral approach. To this purpose, we rewrite Eq. ()22j) with the help of the well 
known trigonometric formula: 

cos((^(t, s) - ip{t, u)) = ^ (e''^(t.^)e-^^(*'") + ^-'^it^s) ^Mt,s)^^ (23) 
Introducing new complex variables 

s) = r (iMe'^(*'") + /le'^i (24) 
Jo 

s) = r dMe-*^(*'") + ke-'^' (25) 
Jo 

(26) 

one finds after simple calculations: 

^'^^^ = 2L Jo "^'^t dT- (2^) 

Now the expression of the energy has become simpler, but the new fields $(s,t) and $(s,t) 
are not independent and have a complicated dependence on the true degree of freedom ip{t, s) 
as shown by Eqs. and (|^. For this reason, it is preferable to assume that $(t, s) and 
^{t, s) are independent complex fields subjected to the constraint jl7|: 

dsHt,s)d.Mt,s) = l (28) 



It is easy to check that Eqs. ()24|1 and ()25|) provide exactly the solution of the constraint 
(I28|) . At this point the energy Ti.o{Lp) may be rewritten in the simple form: 

where $, $ are treated as independent complex degrees of freedom subjected to the condition 
(I28|) . This constraint can be imposed for instance by means of a Lagrange multiplier. 
Instead of the complex coordinates $ one may also exploit real coordinates x,y: 

rs I _ 

x(t,s)= / ducos(fi{t,u) + licosipi = -{^(t, s) + ^{t, s)) 
Jo 2 

y{t,s)= / dusmLp{t,u) + hsiuLpi = —{^{t, s) - ^{t, s)) (30) 
Jo 2i 

It is easy to see that x(t, s) and y(t, s) correspond in the continuous case to the discrete 

coordinates x„ and y„ of the point P„ given by Eq. Q. Written in terms of the real variables 

x{t, s) and y{t, s) the functional Ho{^, $) and the constraint (^Hj) become respectively: 

M 

'^o{x,y) = ^J^ dsix\t,s) + y\t,s)) (31) 
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and 

{dsx{t,s)f + {dMt,s)y = l (32) 
In the discrete case this constraint corresponds to the condition: 

{xi - Xi.i)'^ + (yi - yi^iY = i = 2,...,N (33) 

In the future it will be convenient to exploit the vector notation: 

K{t,s) = {x{t,s),y{t,s)) (34) 

and 

M 

Hoix, y) =no{'R) = — dsR\t, s) (35) 

ZL Jo 

If s = 0, it is clear from Eq. (jHU)) that the point R(t, 0) is fixed for every time t at the 
location: 

R(t, 0) = (Zi cos (fii, li sin (pi) (36) 

This is in agreement with the assumption of Eq. ^2\\ji. where it is supposed that li and ipi 
are constant in time. 

III. THE CLASSICAL EQUATIONS OF MOTION 

We wish to give a "stringy" interpretation of the chain dynamics, in which the chain moves 
during a time interval [tj, t/] from an initial conformation $j(s) to a final conformation $/(s). 
During this motion, the chain spans a two dimensional portion of the plane {x,y), whose 
points are described by the "complex coordinates" $(t, s), $(t, s). The chain "world-sheet" 
is delimited by the range of the variables t and s: ti < t < tf and < s < L. The situation 
is depicted in Fig. |21 Following this strategy, we define the chain action: 

A-^r^r"'^^ (37) 

2LJu Jo dt dt ^ ' 

In the rest of this Section we study the classical equations of motion corresponding to the 
action Ao'- 

~^F~ - - ° ^^^^ 

and the constraint (l2~ 



ds 



1 (39) 
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The solutions of Eqs. ()38|) are of the form: 

^,iit, s) = ais) + this) s) = ais) + this) (40) 

a and b are complex functions of s and a, b are their complex conjugates. If we put 

, , tf^i(s) — ti^f(s) , 1 ,^ , . :^ / XX /,,X 

= ^ \^ ; ' Ks) = -i^fis) - Hs)) (41) 

tf — ti tf — ti 

then Eq. ()4U|) represents the evolution of a chain which during the time tf — ti passes from 

an initial conformation $j(s) to a final conformation $/(s) for < s < L. The constraint 

fl28|l requires additionally that b and b are constants independent of s and that: 

dsais) = e'^^'^ dsa{s) = e-'^^'^ (42) 

where (^(s) describes the angles of a given static conformation of the chain. Going to the 
real coordinates s) and y{t, s) of Eq. (jHUjl . this implies that: 

Xc«(t, s) = Xc/(s) = / (i'UCOS(^(-u) + /i C0S(y9i (43) 

Jo 

Vd {t,s) = yci{s) = du sin i^{u) + /i sin (/^i (44) 
Jo 

We note that in the above equation we have put 6 = 6 = 0, so that there is no dependence on 
the time. This is required by condition ()36p . which demands that the beginning of the chain 
is fixed at the point (/i cos<y9i, li simpi). In other words, in the absence of interactions, the 
conformation of the chain does not change in time. This is not surprising. In fact, we note 
that, in the passage from the discrete kinetic energy Kjf^^ to its continuous counterpart, the 
term Yln=iYlk=i^n-k+i^^n-k+i which was present in Kjf^^ disappeared after performing 
the limit a — > 0. This fact, together with the conditions (PT|) . which fix one of the ends of 
the chain, make the classical dynamics of the ideal chain trivial. 

As anticipated in the previous Section, it is now easy to add the interactions. For example, 
let us suppose that the segments of the chain are immersed in an external potential V"ext(i') 
and that there are also internal interactions associated to a two-body potential Vi„t(ri,r2). 
In this case, Eq. ()27|) generalizes to: 

A = Ao + Ae^t + Ant (45) 

where has been defined in Eq. (|?fj). while 

dt rdsK..(R(t,s)) (46) 
Jti Jo 



rtf 

Aext 
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and 

Aint= r dti [ dSi r dt2 I rfS2Kmt(R-(^l,Sl),R(t2,S2)) (47) 

JU Jo Jti Jo 
Alternatively, if one wishes to express Aext and Amt as functionals of $(t, s) and <i>(t, s) it 

is possible to exploit Eqs. (jHSI). 

IV. DYNAMICS OF A CHAIN IMMERSED IN A THERMAL BATH 

At this point we are ready to study the dynamics of a chain fluctuating in a solution at 
fixed temperature T. Once the action of the system is known, it is possible to introduce 
the dynamics using a path integral approach. To this purpose, we consider the following 
probability distribution: 

= /*(*,.)^.,(.) mt,s)V^t,s)exp (-^) mHt,s)\' - 1) (48) 

where A is the action of Eqs. P3|) - (jTrj) . The boundary conditions for the complex conjugated 
field $ do not appear in Eq. ()48|). because they are fixed by the boundary conditions of 
the field $ due the constraint ()28|1 . which is imposed in Eq. ()48|1 by means of the Dirac 
5— function. Let us also note that due to that constraint the action A may be written in a 
form which resembles more that of a two dimensional field theory: 

M f^f 

^=W7 dt ds\V^{t,s)\^ + Aext + Aint (49) 

2L Jti Jo 

with V = {dt, dg). The actions Aext and Aint are useful in order to describe the interactions 
of the chain with itself and with the surrounding environment. The distribution "^^"^ gives 
the probability that a chain of length L starting from an initial conformation $j(s) at time 
t = ti arrives to a final conformation $/(s) at the instant t = tf after fluctuating in a thermal 
bath held at constant temperature T. To understand how this probability is computed by 
means of the path integral appearing in the right hand side of Eq. PHj) it is helpful to look 
at Fig.nQ- 

During its motion in the time interval [ti, tf] the chain spans a surface S. In 
Eq. (PH|) it is performed a sum over all possible surfaces of this type. If the temperature T 
is zero, the chain moves according to the classical equations of motion. In that case, the 
conformations of minimal energy are favored and the energy is conserved along the motion 
at each instant. If however the temperature is different from zero, the energy does not 
need to be conserved in time, because of the thermal fluctuations. The conformations for 
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FIG. 2: In its motion from the initial conformation ^j(s) to the final conformation $/(s) the chain 
spans a two dimensional surface in the space. 



which the energy is not minimal are suppressed in the exponential exp 



A 



appearmg 



in Eq. (j^8|) . but may still give a relevant contribution to the whole path integral if their 
number is overwhelming with respect to the conformations of minimal energy. 

Let us now discuss the boundary conditions which appear in the path integrals of Eq. ()48j] . 
Since the interactions are not relevant in the present context, we will consider just ideal 
chains. It is also convenient to reformulate Eq. ()48|) in terms of the real variables x{t, s) and 
y{t,s): 



2d 



R(ti,s)=Ri(s) 



Vx{t, s)Vy{t, s)e "bt § ((dsxf + {d.y) 



where 



M fif rL 



2L Jti Jo 



Apart from the boundary conditions: 



R{ti,s) 



Rfis) 



(50) 
(51) 

(52) 
(53) 



and 



R(t, 0) = (/i cos ipi, li sin Lpi) (54) 

one could also add the requirement that the chain forms a closed loop. This further condition 
is implemented by the constraint: 



R(t,0) = R(t,L) 



(55) 
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Eq. ()55|) constrains the two ends of the chain, occurring when s = and s = L, to coincide 
with the fixed point (/i cos (^i, /i sin (/^i) at any instant t. It is cumbersome to keep track 
of all these boundary conditions within the path integral, so that it becomes preferable to 
expand the fields R(t, s) around a classical background: 



R(t,s) = Re/(t,s) + R,(t,s) 



(56) 



We first concentrate ourselves on the background. Due to the fact that there we assumed 
that the chain is ideal, the classical fields Rc;(t, s) can be chosen among those which obey 
the free equations of motion associated to the action Aq of Eq. (|^: 



Relit, s) = 



(57) 



The above equation should be solved in such a way that its solutions satisfy the bound- 
ary conditions ()52|) - (jHH), but not the constraint ()28|) . because this is taken into account 
separately by the 6 functions appearing Eq. (jKIH) . Explicitly, one may write the general 
expression of Rc/(t, s) as follows: 



(58) 



tf — ti tf — ti 

The boundary condition (jH^ is fulfilled by requiring that Ri(0) = R/(0) 
(/i cos v^i, /i sin (/Ji). If one wishes to add closed loop condition (fS^. one has to ask ad 
ditionally that: 



R^(0) = Rr{L) = Ro(0) =Ro(L) = {h cos(^i,/isin(^i) 



(59) 



Let's now discuss the component Rq{t,s) appearing in Eq. (jSUj) . It describes the statistical 
fiuctuations around the classical background Rc/(t, s). Since the non-trivial boundary con- 
ditions are already taken into account by the background, Rq{t, s) has a trivial behavior at 
the boundary: 

R,(ti,s) = R,(t/,s) = (60) 

After having split the fields R(t, s) as in Eq. (jSUj) . the probability distribution (I^Uj) be- 
comes: 

2d „ \ -4,0 (Rd) 



exp 



X j VRqit, s) exp 



A(Rg) 



6 {{d^Rli) + (dsRl) + 2dsRei ■ dsRq - 1 



(61) 
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where Aoijlci) contains just the classical field conformations: 

M rL 



M 



while the statistical fiuctuations are in A()(Il„ 



AoCRq 



M 
2L 



dt 



dsRi 



(62) 



(63) 



After introducing a Lagrange multiplier X{t, s), it is possible to write in Eq. ()61|) as follows: 



2d 



exp 



Ao(Rci) 



X J V'Rg{t, s) exp 



^o(R-g 



X exp 



dt / ds 







(64) 



-A(t, s) ((9,R,0' + {dsK,y + 2dsKci ■ dsKg - l] 

After the field splitting of Eq. (jSUj) . one has to deal only with the statistical fiuctuations 
Rq, whose boundary conditions are trivial. In this way it becomes in principle possible to 
eliminate the Rg's from the path integral ()64|1 by performing a gaussian integration. Finally, 
let us discuss when it is convenient to split the fields as in Eq. ()56|1 . We know already from 
Section Uni that the constraint forces the identity R/(s) — Ri{s) = in Eq. (jKHjl as 
in the case of the classical solutions ()43|) and For this reason, if the final and initial 

conformations R/(s) and Rj(s) are chosen to be very different, it is licit to expect that 
the statistical fluctuations Rg{t,s) will not be small, since they are needed to restore the 
constraint p8|) . Thus, from a perturbative point of view, it is possible to consider the 
fluctuations Rq{t, s) as small perturbations only if R/(s) ~ Ri(s), i. e. 

|R/(s) -R,(s)| < e (65) 

with e being a small constant parameter. 



V. THE THREE DIMENSIONAL CASE 



In this Section we consider a chain of segments PiPi-i of flxed lengths k, i = 2, . . . , N, 
in three dimensions. Using spherical coordinates, the positions of the end points of the 
segments Pi{t) = {xi{t),yi{t), Zi{t)) are given by: 

Xn{t) = h cos ipi{t) sin 9i{t) 

Vnit) = Elilismipi{t)smeiit) n = l,...,N (66) 

Znit) = Etik cos 9i{t) 
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For simplicity, the kinetic energy K^fg^ of the system will be computed in the particular case 
in which the chain is attached at the origin of the coordinates, i. e. 

Pi = (0,0,0) l^ = i^ = Q (67) 

After a long but straightforward calculation, the result is: 

^disc = X! X! l^^n-k+lVn-k+l sin 9n-k+l + X! X! ~l^^n-k+l^n-k+l 
n=l k=l ^ n=l A:=l ^ 

A' 71—1 n— fc 

+ XI XI X! '^nkL^k+l [pi'pn-k+l sin 6'j sin 6'„_A:+1 cos (V5n-fc+l - V'i) 
n=l k=l 1=2 

+ COS Bi sin ^„_A:+1 sin (</?» - ^n-k+l) 

+ ipiOn-k+i sin 6'i COS 6'„_a:+i sin (v?„-fc+i - <^i) 

+ Ofin-k+i (cos 6'i COS 6'„_fc+i COS (v5n-fe+i " V^i) + siu 6*^ sin dn-k+i) (68) 

It is now possible to pass to the limit of a continuous chain a — and — > +oo. Making 
the same assumptions of uniform length and mass distributions as in Eqs. fjl9p and ()20|). we 
get after some algebra: 
M fL 



j^3d _ 

I 



/ ds{L — s) / dv [(p{t,v)(p{t, s) sin 9{t,v) sin 9{t^ s) cos {(p{t, s) — (p{t,v)) 
Jo Jo 

+ 9(t, v)ip{t, s) cos 9(t, v) sin 9{t, s) sin {(p{t, v) — ip(t, s)) 
+ 9(t, s)ip{t, v) sin 9{t, v) cos 9{t, s) sin {(p{t, s) — ip{t, v)) 

+ 9{t, v)9{t, s) (cos 9{t, v) cos 9{t, s) cos ((^(t, s) - (p{t, v)) + sin 9{t, v) sin 9{t, s))] (69) 



This is the three dimensional analog of Eq. 

The kinetic energy of (j69|) is complicated, but from the lesson of the two dimensional 
case we know how to simplify it. First of all, we note that the kinetic energy of the discrete 
chain may be written in terms of the cartesian coordinates (j66p as follows: 

KlL = T.'^i^l + yl + ^l) (70) 

n=2 ^ 

where a;„, y„ and Zn have been defined in Eq. ()66|1 . The sum over n starts from 2 because 
one end of the chain coincides with the origin of the axes, so that li = 0. Of course, due to 
the condition that each segment has a fixed length Eq. (f7n|) must be completed by the 
following constraints: 

{Xn - Xn-lY + iVn - Vn-lf + {Zn - Zn-lf = li H = 2, . . . , N (71) 
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At this point we have two choices. Either we keep the kinetic energy in the simple form 
of Eq. (f7n|) at the price of having to deal with the constraints ()71|) . or we solve those 
constraints using spherical coordinates li,6i,(j}i, as it has been done in Eq. (jMjl . In the latter 
case, we have the complicated expression of the kinetic energy of Eq. and there is the 
problem of defining a path integration over spherical coordinates. In the continuous limit, 
the situation does not change substantially. After performing the continuous limit following 
the prescriptions of Section |Tll the kinetic energy of Eq. (fTOj) and the constraints (fTTjl are 
respectively replaced by: 



2L Jo 



idtx{t, s)f + (dtvit, s)f + (dtzit, s)f (72) 



and 

{dMt, s))' + {OMt, s)f + {dsz{t, s)y = l (73) 
The constraint (f73|) can be eliminated by introducing spherical coordinates 6{t, s), ip{t, s): 

x{t,s) = / du cos (p{t,u) sin 9 {t,u) (74) 

<^ 

y{t,s) = / du simp {t,u) sin 6 {t,u) (75) 

<^ 

z{t,s) = / ducos6{t,u) (76) 



If one makes the coordinate substitutions of Eqs. (fTIHTm) in the kinetic energy (ff^ and 
applies the formula ()17|). one arrives exactly at the expression of the kinetic energy (jU^ . 
Thus, Eq. flB^ and Eq. (f7^ together with the constraint (f7H|) are equivalent. 

To construct of the probability distribution in three dimensions we choose the ap- 
proach in which the the coordinates x{t, s), y{t, s) and z{t, s) are independent and the right 
number of degrees of freedom is restored by the condition ()73|) . The result is similar to that 
of the two dimensional case: 

= /h,,.)^k,,) ^nt, s) exp |-^| sm^r - 1) (77) 



where 



•Af = ^ dt ds\{dtx{t,s)f + {dty{t,s)f + {dtz{t,s)f\=— dt dsB? (78) 

IL Jti Jo '- -' ZL Jti Jo 

The above expression of the probability distribution has been obtained assuming that one 
end of the chain is attached at the origin according to the condition of Eq. ()67|). However, 
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using translational invariance, it is easy to check that Eqs. (fTTjl and (fTSjl remain vahd also 
if the chain has one of its ends fixed at any other given point Hfi^ed, so that to complete 
the set of boundary conditions of Eq. (|77jl we may add the following requirement: 

R(t,0) = R,(0) = R/,,ed(0) (79) 

R'/ia;ed(0) being a fixed point. 



VI. CHAINS WITH CONSTANT ANGLE OF BENDING 

The approach presented above in order to treat the dynamics of random chains has some 
interesting variants which we would like to discuss in this Section. To this purpose, we 
choose the formulation in which the positions of the ends of the segments composing the 
chain are given in cartesian coordinates. As we have already seen, in this way the expression 
of the kinetic energy K^f^^ is simply: 

K!L = j:^iil + yl + il) (80) 

n=l ^ 

However, one has to take into account also the constraints: 

(Xn - Xn-lY + iVn - Vn-lY + - Zn-lY = ll U = 2, . . . , N (81) 

We assume as before that all segments have the same fixed length = a, but additionally 
we require that: 

(z„ - 2„„i)2 = b^<a^ (82) 

This implies that the projection of each segment onto the 2;— axis has length ±6, so that the 
segments are bound to form with the 2;— axis the fixed angles ai = a or a2 = {ti ~ a) defined 
by the relations: 

b b 

costti = H — cosa2 = — (.83) 

a a 

Clearly, in both cases the constraints (jHT|) and (jH^ may be rewritten as follows: 

{Xn - Xn-lf ^ iVn ' Vn-lf ^ _ ^ n = 2, . . . , iV (84) 



cos^ a 

We suppose now that only the angle ai is allowed, so that the chain cannot make turns in 
the z direction. An example of a conformation of a chain satisfying these assumptions is 



z 
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FIG. 3: Example of motion of a chain whose segments are constrained to form a fixed angle a with 
the z— axis. In the figure a = 30° 

given in FiglHl The constraints ()84j] are solved by choosing spherical coordinates, in which 
however the angles On formed by the segments with the 2;— axis are always equal to a: 

n 

^nif) = ^ /j COS sincK (85) 

i=l 
n 

yn{t) = ^ h sin ipi{t) sin a (86) 
1=1 

n 

^n{t) = ^ /i COS a = n cosa (87) 
1=1 

As we see from the above equation, each segment is left only with the freedom of rotations 
around the direction, corresponding to the angles ifi{t). Moreover, the total length of the 
chain is always L = Na, but now also the total height h of the trajectory along the 2;— axis 
is fixed: 

h = Nb (88) 

At this point, we pass to the continuous limit, this time taking as parameter describing the 
trajectory of the chain the variable z instead of the arc- length s. Due to the last of Eqs. (jHTj) . 
the 2;— components of the velocities are always zero: 



i„(t) = 



(89) 
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As a consequence, we are left with a two dimensional problem, which may be treated in 
exactly the same way as we treated the two dimensional chain of Section The only 
difference is that Eqs. Q should be replaced by Eqs. (jHHjl and (jHHjl and the constraints have 
a slightly different form. Following the same procedure presented in Section |n] we find after 
a few calculation the expression of the kinetic energy: 

= tan^a dz dzi / dz2Pm{t, z)pi{t, z - Zi)pi{z2) 

Jo JQ Jo 

X ip{t,Z - Zi)(p{t,Z2)cOs{(p{t,Z - Zi) - (p{t,Z2)) (90) 

and of the constraint ()84|1 : 

{d,xf + {d,yf = tan^ a (91) 

Apart from the appearance of the factor tan^ a and the choice of the height z instead of the 
arc-length s, Eqs. ^ and ^ are identical to Eqs. jHI) and (HH) in the limit l{t) = i{t) = 0. 
It is now not difficult to show that the probability distribution is given by: 

^l' = J Vx{t,z)Vy{t,z)exp^^-^^6{{d,xf + {d^y f-te.n'a) (92) 



X + y 



(93) 



where 

Aoa = tan^ a [ dt [ > 
Jti Jo 

At this point we discuss briefly the case in which both angles ir — a and a are allowed. 
In this situation, the trajectory of the chain may have turns. An example of motion of this 
kind is given in FiglU The constraints ()8H1 and ()82|1 remain unchanged, but the coordinate 
z cannot be chosen as a valid parameter of the trajectory of the chain and one has to 
come back to the arc-length s. The most serious problem occurs due to the fact that the 
variables Zn{t) are not continuous functions of the time, since each Zn{t) is allowed to jump 
discretely between the two discrete values +b and —b, corresponding to the angles a and 
TV — a respectively. It is therefore difficult to define the components i„ of the velocities of 
the ends of the segments and thus their contribution to the kinetic energy. Let us note that 
this problem affects only the z degrees of freedom. The degrees of freedom Xn{t) and ?/„(t) 
of the chain remain continuous functions of t despite the jumps of the z„'s. This fact can 
be easily verified looking at the definition of x„(t) and ?/„(t) in Eqs. (jH3j) and (jHUj) . Since 
sin(7r — a) = sin a, both the x„(t)'s and yn(t)'s are not affected by the jumps of the angle 
a < — > n — a. The situation simplifies only if the chain has no interactions in which the z 
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FIG. 4: Example of motion of a chain whose segments are constrained to form fixed angles a. or 
vr — a with the z— axis. In the figure a = 30°. Turning points are emphasized by means of black 
points. 



variable is involved. The reason is that in the kinetic energy and in the constraints given 
by Eqs. (jHUHH^ the degrees of freedom connected to the motion along the z— directions 
are decoupled from the other degrees of freedom and may be ignored. As a consequence, 
in absence of z— dependent interactions, the difficulties related to the motion along the 
2;— direction disappear and once again the problem reduces to that the two dimensional 
chain treated in Section |nl Since the constraints are always those of Eqs. (jHDHH21) one may 
proceed as in the case of fixed angle a. As a result, one finds that the final probability 
distribution is of the form: 



C J Vx{t, s)Vy{t, s) exp 



where 



0,Q,7r— a 



0,ci,7r— a 



sm^ a [ dt J ds 



X + y 



(94) 



(95) 



and C is a constant containing the result of the integration over the decoupled z degrees of 
freedom. Let us note in Eqs. ()94|) and ()95|) the appearance of the factor sin^ a in the action 
instead of tan^ a and the replacement of z with the arc-length s as the parameter of the 
trajectory of the chain. 
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VII. CONCLUSIONS 

The dynamics of a random chain has been discussed both from the classical and statistical 
point of view. The classical energy of the chain has been derived in two and three dimensions. 
We have mainly concentrated ourselves on the computation of the kinetic energy because 
the addition of the contribution of the interactions to the total energy is straightforward, see 
Eqs. fl46|) and (07j). The kinetic energy on the contrary has a complicated and non-Markoffian 
expression as shown by Eq. (jH)) in two-dimensions and by Eq. ()68|) in three dimensions. After 
passing to the continuous limit discussed in Section II, one term disappears from the kinetic 
energy. This is the reason for which in the continuous case the classical equations of motion 
have simple solutions of the form given in Eqs. ()43j) - P^ . Let us note that the passage to 
the continuous limit is straightforward and does not need mathematical subtleties as the 
analogous limit by which the Edwards model is obtained in the statistical mechanics of 
random chains. 

The stringy approach to the dynamics of random chains is formulated in Section IV. One 
would naively expect that the final model is related with non-relativistic string theories, 
such as those derived in Refs. but this is not the case. In two and three dimensions the 
probability distributions and \Ef^'^ are respectively given in Eqs. ()48|) and ()77 p - (f7Hj) . As it 
is possible to see from these equations, the path integral sums which provide the expressions 
of \E'^'^ and \E''^'^ have the form of a 0{n) nonlinear sigma model on a two dimensional world- 
sheet, where n = 2, 3 depending on the dimensionality of the physical space in which the 
chains fluctuate. The difference with respect to the nonlinear sigma model is that here 
the derivatives of the fields with respect to the arc-length s are subjected to the condition 
(9sR)^ = 1. In standard nonlinear sigma models it is instead the modulus of the fields 
themselves to be constrained. The presence of a complicated constraint and the absence 
of a small parameter which could be used to start a perturbative expansion complicate the 
computation of the probability functions. For this reason, in Section IV it has been proposed 
the field splitting This splitting provides a convenient way to deal with the boundary 

conditions satisfied by the fields R(t, s) and also allows a perturbative treatment provided 
the initial and final conformations Ri(s) and R/(s) do not differ very much in the sense 
of Eq. (jnH). Alternatively, the statistical fluctuations Rq(t,s) may be eliminated from the 
partition function of Eq. ()64|1 with a Gaussian integration |l^ , but then one ends up with 
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a complicated two dimensional field theory in which the Lagrange multipliers X{t, s) are 
coupled together with the classical background Rci(t, s). Another possibility in order to 
simplify the model consists in relaxing the constraint (S^R)^ = 1, imposing for instance the 
weaker condition lu lo rfs(9sR)^ = 1. This constraint requires that the average length 
of the chain in the interval of time ti < t < tf is L. In this case the Lagrange multiplier A 
does no longer depend on t and s and calculations become easier. It would also be interesting 
to study the probability functions \I'^'^ and in the limit tf — > +oo and tj — > — oo or 
to sum them over all possible values of the initial and final conformations Ri(s), R/(s). 

Finally, chains with fixed angles have been discussed in Section VI. Our approach is valid 
only if the chain has no turning points. If there are turning points the kinetic energy is 
not well defined, because the variable z{t, s) is no longer a continuous function and thus its 
time derivative becomes a distribution. One way for adding to our treatment turning points 
as those of Fig. |3] is to replace the variable z with a stochastic variable which is allowed 
to take only discrete values. Another way is to look at turning points as points in which 
the chain bounces against an invisible obstacle. A field theory describing a one-dimensional 
chain with such kind of non-holonomic constraints has been already derived in Refs. jl^ . 
The problem of turning points is currently work in progress, as well as the possibility of 
imposing topological constraints on the trajectory of the chain. 
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